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Purpose: Systemic lupus erythematosus (SLE) is an autoimmune connective tissue disease that can affect the cen- 
tral nervous system. Neuropsychiatric symptoms are found in 25-70% of patients. Using diffusion tensor imaging 
(DTI) various studies have reported changes in white matter integrity in SLE patients with neuropsychiatric 
symptoms (NPSLE patients). The purpose of this study was to investigate, if regional changes in white matter in- 
tegrity can also be detected in SLE patients without neuropsychiatric symptoms (non-NPSLE patients). 
Methods: Applying DTI and tract based spatial statistics (TBSS) we investigated 19 NPSLE patients, 19 non-NPSLE 
and 18 healthy controls. Groups were matched for age and sex. Image pre-processing was performed using FSL, 
following the TBSS pipeline (eddy current correction, estimation of fractional anisotropy (FA), normalization, 
skeletonization of the group mean FA image). A general linear model with threshold-free cluster enhancement 
was used to assess significant differences between the three groups. 

Results: Statistical analyses revealed several regions of decreased prefrontal white matter integrity (decreased 
FA) in both groups of SLE patients. The changes found in the non-NPSLE patients (as compared to healthy con- 
trols) overlapped with those in the NPSLE patients, but were not as pronounced. 

Conclusions: Our data suggest that changes in regional white matter integrity, in terms of a decrease in FA, are 
present not only in NPSLE patients, but also in non-NPSLE patients, though to a lesser degree. We also demon- 
strate that the way statistical maps are corrected for multiple comparisons has a profound influence on whether 
alterations in white matter integrity in non-NPSLE patients are deemed significant. 

© 2014 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license 

(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Systemic lupus erythematosus (SLE) is a clironic inflammatory, 
immune-mediated disease, whicii affects 0.05-0.1% of tlie general pop- 
ulation (Somers et al., 2014; Bernatsl<y et al., 2007). Neuropsychiatric 
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dMRI, diffusion MRl; DTI, diffusion tensor imaging: FA, fractional anisotropy: MRl, magnet- 
ic resonance imaging; SLE, systemic lupus erythematosus: SLEDAl, Systemic Lupus 
Erythematosus Disease Activity Index; SLICC, systemic lupus erythematosus 
International Collaborating Clinics: SD, standard deviation; SVM, support vector machine; 
TBSS, tract based spatial statistics: TFCE, threshold free cluster enhancement: VBM, voxel 
based morphometry. 
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symptoms have been reported to occur in 25-70% of patients with 
SLE (Brey et al., 2002; Monastero et al., 2001) and are associated 
with increased morbidity and mortality (Sibley et al., 1992; 
Gladman et al., 1996). Magnetic resonance imaging (MRl) has be- 
come part of the routine clinical work up in SLE patients presenting 
with neuropsychiatric symptoms (NPSLE patients) (Appenzeller et al., 
2008a,b; Steup-Beekman et al., 2013), helping to diagnose or exclude ce- 
rebral pathologies such as edema, hemorrhage, central thrombosis, or 
stroke. 

Using quantitative analysis methods, changes in brain chemistry 
(Sundgren et al., 2005; Appenzeller et al., 2007a), regional cerebral 
blood flow (Shen et al., 1999; Oda et al., 2005; Wang et al., 2012) and 
cerebral glucose utilization (Otte et al., 1997) have been demonstrated 
in NPSLE patients. Structural brain imaging methods, such as voxel 
based morphometry (VBM), cortical thickness analysis and diffusion 
tensor imaging (DTI), have also been applied to investigate differences 
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in regional gray and wliite matter in SLE patients and healtliy controls 
(HCs) (Appenzeller et al., 2007a; Jung et al., 2010; Emmer et al., 2010; 
Welsh et al., 2007). Regions of altered brain structure, in terms of de- 
creases in regional gray matter volume, cortical thickness or fractional 
anisotropy (FA) have been demonstrated in NPSLE patients, but less fre- 
quently in non-NPSLE patients (Appenzeller et al., 2007a; Jung et al., 
2010; Jung et al., 2010). 

The distinction between NPSLE and non-NPSLE is clinical, based on 
the American College of Rheumatology (ACR) case definitions. In clinical 
practice brain imaging is often not part of the routine work up of non- 
NPSLE, possibly leading to an underestimation of CNS disease. Some au- 
thors, using MR imaging, have found evidence for small emboli, cerebral 
infarctions, hemosiderin deposits, and dilated perivascular spaces in up 
to 60% of non-NPSLE patients, strongly suggesting cerebral involvement 
in at least some of these patients (Cotton et al., 2004). Along these lines 
spectroscopy studies have reported decreased NAA levels, indicative of 
axonal dysfunction/damage, in SLE patients with active disease, inde- 
pendently of CNS manifestation (Sundgren et al., 2005; Appenzeller 
et al., 2005; Cagnoli et al., 2013). On the other hand some of the 
studies using advanced statistical imaging methods, such as VBM and 
DTI, have detected reduced gray or white matter integrity in NPSLE pa- 
tients only, but not in non-NPSLE patients, suggesting that non-NPSLE 
patients with respect to gray and white matter structure are more sim- 
ilar to HCs than to NPSLE patients (Appenzeller et al., 2007b; Jung et al., 
201 Oa,b). For example, Jung et al. investigating NPSLE patients and non- 
NPSLE patients could identify regions of decreased white matter integ- 
rity (decreased FA values) in NPSLE patients only (Jung et al., 2010b). 
However, in statistical neuroimaging other factors than "true" differ- 
ences (of outcome parameters, e.g. gray matter density/volume, FA 
values, etc.), such as group size, acquisition parameters, preprocessing 
steps and the way multiple comparisons are accounted for play a role 
for whether differences are considered and reported as significant or 
not. 

The aim of this study was to further investigate white matter integ- 
rity in NPSLE and non-NPSLE in a cohort of 19 NPSLE patients, 19 non- 
NPSLE patients and 18 HCs. We used FSL and the TBSS pipeline (tract 
based spatial statistics) to preprocess and statistically analyze our data 
performing both univariate and multivariate statistical tests. We hy- 
pothesized that the pattern of altered white matter integrity in non- 
NPSLE patients does not substantially differ from that seen in NPSLE pa- 
tients. We also sought to specifically investigate the impact of alterna- 
tive ways to correct for multiple comparisons that are frequently used 
in statistical brain imaging. We hypothesized that the way statistical 
maps are created has a profound influence on whether non-NPSLE pa- 
tients have significant alterations in white matter integrity. 

2. Material and methods 

2.1. Subjects 

The study comprised 1 9 patients with acute NPSLE, defined as one or 
more neuropsychiatric manifestations within 2 weeks prior to inclusion 
in the study, and of 1 9 patients with SLE, without current or past histoiy 
of neuropsychiatric symptoms (non-NPSLE), as well as of 18 HCs. In a 
recently published VBM study, we had reported differences in regional 
gray matter volume between NPSLE patients, non-NPSLE patients and 
HCs (Cagnoli et al., 2012). Some of these patients/subjects had also un- 
dergone DTI and took part in current study. There were 19 female pa- 
tients (out of 19) in the NPSLE group, 18 female patients (out of 19) in 
the non-NPSLE group, and 16 females (out of 18) in the HC group. The 
gender ratios did not differ significantly between groups (Fisher's 
exact test, p = 0.32). 

The classification of NPSLE was based on the 19 ACR case definitions 
(American College of Rheumatology, 1999). In order to participate in 
the study, individuals needed to fulfill 4 or more of the ACR classification 
criteria (Tan et al., 1982) and, in order to qualify as an NPSLE patient. 



present within 2 weeks of the initial neuropsychiatric event. Patients 
in the non-NPSLE group were required to have no prior history of neu- 
ropsychiatric manifestations regardless of attribution. All participants 
were recruited from the Department of Internal Medicine, Division of 
Rheumatology, and the Department of Emergency Medicine at the Uni- 
versity of Michigan. Patients with previous/known drug abuse, alcohol 
abuse, diabetes, stroke and/or renal insufficiency were excluded. Groups 
were matched for age and sex. For details on patient groups and clinical 
manifestation see Tables 1 and 2. The study had been approved of by the 
IRB of the University of Michigan (HUM00000714) and informed con- 
sent had been obtained from the patients prior to study inclusion. 

All patients and HCs were submitted to a standardized clinical and 
neurological examination including medical history (age, gender, eth- 
nicity, age at disease onset, disease duration, etc.), physical examination 
standard laboratoiy assessment and a Mini Mental State examination 
(Folstein et al., 1975). SLE and NPSLE patients unden/vent required 
laboratory tests to determine disease activity by the Systemic Lupus 
Erythematosus Disease Activity Index (SLEDAI) (Bombardier et al., 
1992) as well as antiphospholipid antibody status. SLE was considered 
active with a SLEDAI score >6. We chose this higher SLEDAI level to 
minimize the inclusion of patients who have persistently active serol- 
ogies (anti-DNA antibodies, low complements and /or leukopenia) 
without clinical evidence of disease activity. 

Disease duration was defined as the time between the diagnosis of 
SLE and the day of the MRl. CNS manifestations were divided into cen- 
tral and peripheral following ACR case definitions. Cumulative damage 
in both patient groups was assessed using the Systemic Lupus Interna- 
tional Collaborating Clinics/American College of Rheumatology Damage 
Index (SLICC) (Gladman et al., 1996) at the time of MRl acquisition. 
With respect to age at disease onset, disease duration, SLEDAI, SLICC 
and Mini-mental scores, non-NPSLE patients and NPSLE patients were 
compared using the Mann- Whitney U test. Fisher's exact test was ap- 
plied to compare gender and ethnic distribution in the three groups 
(non-NPSLE, NPSLE and HCs). For details see Table 1. 

Patients and HCs with MRl pathologies that affected white matter in- 
tegrity, such as cerebral infarcts or confluent white matter lesions, were 
excluded from the study. Other incidental findings that did not affect 
white matter did not lead to an exclusion (for a list of incidental findings 
see Section 3.2). Small, unspecific white matter lesions were found in 
both patients and HCs ( see Results), they were not considered patholog- 
ical per se, as they are found in up to 45% of healthy people (Fazekas 
et al., 1988; Schmidt et al., 1992), and did not lead to an exclusion of 
the subject. 

2.2. Data acquisition, pre-processing and analysis 

MR imaging was performed at a 3 T MR scanner (Achieva, Philips, 
Best, Netherlands). To rule out gross abnormalities such as acute or 
old infarcts, hemorrhage, focal lesions or pathological contrast enhance- 
ment and for any additional lesions not related to SLE, T2-w, FLAIR, Tl- 
weighted images pre and post contrast administration were acquired in 
the same session and evaluated by an experienced neuroradiologist. 
Lesion burden in form of white matter T2/FLA1R hyper-intensive lesions 
was defined as mild (1-5 small white matter lesions), moderate (5-10 
white matter lesions), or severe (>10 white matter lesions). Atrophy 
was visually evaluated as none, mild, or severe. 

dMRl images (15 gradient directions atb = 800 s/mm^, plus one 
b = 0 image) of 19 NPSLE patients, 19 non-NPSLE and 18 healthy con- 
trols (HCs) were acquired. Image pre-processing and statistical analyses 
were performed using FSL 5, following the standard TBSS pipeline: eddy 
current correction, linear regression of the diffusion tensor model, FA 
estimation and non-linear normalization (FNIRT) were followed by 
skeletonization of the group mean FA image. Thresholding at FA > 0.2 
to restrict analysis to the white matter led to 122,477 mean skeleton 
voxels, to which data from individual subjects was projected. A general 
linear model was used to assess significant differences between the 
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Table 1 

Demographic data and disease severity. 





NPSLE patients 
n = 19 


Non-NPSLE patients 
n = 19 


HC 

n = 18 


Significance 




Range 


Range 


Range 


p value 




mean/SD 


mean/SD 


mean/SD 




Age 


22-67 


21-60 


18-61 


Kruskal-Wallis test 




39.1/11.7 


38.0/11.1 


37.6/13.5 


p = 0.99 


Female patients 


19 


18 


16 


Fisher's exact test 










p = 0.32 


Ethnicity 


1/19 Asian 


2/19 AA 


1/18 Asian 


Fisher's exact test 




5/1 9 AA 


17/19 White 


17/18 White 


p = 0.072 




13/19 White 








Age at disease onset (years) 


12-65 


20-55 


n/a 


Mann-Whitney U test 




33.4/13.0 


33.3/10.8 




p > 0.99 


Disease duration (years) 


0.1-16 


1.0-15 


n/a 


Mann-Whitney U test 




5.7/6.3 


4.7/4.8 




p = 0.95 


SLEDAl 


1-24 


0-8 


n/a 


Mann-Whitney U test 




10.8/5.5 


2.0/2.2 




p < 0.0001 


SLICC 


0-2 


0-3 


n/a 


Mann- Whitney U test 




0 68/0.67 


0.47/0.84 




p = 0.18 


Mini Mental 


22-30 


24-30 


n/a 


Mann-Whitney U test 




27.8/2.6 


28/1.5 




p = 0.60 



AA = African American; HC = healthy contr ol; NPSLE = neuropsychiatric systemic lupus erythematosus; SD = standard deviation; SLE = systemic lupus erythematosus; SLED/U = Systemic 
Lupus Erythematosus Disease Activity Index; SUCC = Systemic Lupus Erythematosus International Collaborating Clinics; SD = standard deviation. 



three groups (all SLE vs. HC; NPSLE vs. HC; non-NPSLE vs. HC); age was 
included as a covariate of no interest. 

To specifically investigate the impact of correcting for multiple com- 
parisons we applied three different methods. First we applied threshold 
free cluster enhancement (TFCE), an alternative approach for control- 
ling family wise error, implemented and recommended in version 5 of 
FSL (Salimi-Khorshidi et al., 2011; Smith and Nichols, 2009). TFCE is a 
voxel based method; it resembles cluster-based analysis in that it boosts 
the confidence in the significance of differences that are observed con- 
sistently over large spatial neighborhoods, but avoids having to select 
an initial cluster-forming threshold (for more details see below). Sec- 
ondly we applied the method used in the paper by Jung et al. to facilitate 
comparison, thresholding the t-statistic images from a General Linear 
Model at f > 3.0, and testing the resulting clusters for statistical signifi- 
cance based on 10,000 Monte Carlo permutations for each group differ- 
ence. We also performed the same procedure applying an initial voxel 
threshold of t > 2.0, to further elucidate the impact of this initial thresh- 
old on cluster statistics. 

We finally performed a multivariate analysis, based on support 
vector machines (SVMs). Specifically, we were interested whether we 
could train a classifier that could differentiate between the three groups 
by identifying patterns of altered FA maps in SLE patients. Even if 
classification accuracy is below the level that would be required for 



Table 2 

Differences in regional fractional anisotropy (TFCE corrected). 



Comparison of groups 


Cluster size 


Coord 


inates 




Smallest p value 








X 


y 


z 




2 


NPSLE < HC 
















7291 


24 


25 


1 


0.008 




Non-NPSLE < HC 












3 




2560 


-27 


37 


2 


0.024 






550 


-17 


0 


38 


0.038 


4 




258 


18 


34 


12 


0.039 




NPSLE < non-NPSLE 


























3. 



HC = healthy control, NPSLE = neuropsychiatric systemic lupus erythematosus; non- 
NPSLE = systemic lupus erythematosus without neuropsychiatric symptoms. TFCE = 
threshold free cluster enhancement. Per-voxel p values, coll ected for mulriple compari- 
sons, were computed using TFCE and permutarion based statistics. Coordinates refer to 
voxels with smallest p values within a cluster. Unlike cluster-based statistics, per-voxel 
testing can lead to many small clusters of significant voxels. For brevity, only the ones larg- 
er than 100 voxels are reported. 



computer-aided detection of the disease, classification rates that are sig- 
nificantly greater than random chance confirm that significant differ- 
ences exist between the successfully classified groups. Besides in cases 
where only one of the two alternative methods for controlling family- 
wise errors indicates a significant difference between two groups, suc- 
cessful classification can serve as an additional confirmation that this ef- 
fect can be attributed to improved sensitivity, rather than to a false 
positive detection. For classification, we used projected FA values at all 
122,477 skeleton voxels as a feature vector. We averaged the FA values 
of each subject in the training dataset over the skeleton and used linear 
regression on the resulting averages to estimate the effect of age. The 
resulting slope and intercept were used to age-correct all feature values. 
In order to reduce the impact of voxels less relevant to the classification 
task, each feature was weighted by its f score, as it would be used in a 
two-sample F test. Classification accuracy was estimated by training a 
linear SVM on the resulting vectors in a leave-one-out manner; care 
was taken to exclude test data from linear regression and feature 
weighting via the F score. Multivariate analysis was performed within 
the shogun toolbox ( http://shogun-toolbox.org/ ). We used permutation 
testing to check if accuracies were significantly higher than random 
chance. To this end, we built up a null distribution by computing 
leave-one-out estimates of the accuracy of classifiers that were trained 
on the same data with permuted labels. Overall our analyses comprised 
the following steps: 

1 The threshold free cluster enhancement (TFCE) method as implemented in FSL 5 which 
is based on pre-processing the t-statisric images using a nonlinear image transforma- 
tion that boosts maxima with support from a spatial neighborhood before performing 
a voxel-based permutation test to correct for multiple comparisons. 
A cluster based method where group FA differences are calculated by thresholding the 
t-statistic images at a cluster-forming threshold of t > 3.0 and controlling for multiple 
comparisons on a cluster level using permutation methods with FSLs Randomize, run- 
ning 10,000 Monte Cario permutations for each of the group differences. 
The same cluster-based method as in 3 with an initial cluster-forming threshold of t > 
2.0. 

A multivariate analysis to train a classifier that could differentiate between the three 
pairs of groups by identifying patterns of altered FA maps. 

3. Results 

3.1. Patients 

Nineteen NPSLE and 19 non-NPSLE patients, as well as 18 HCs partic- 
ipated in the study. The three groups were matched by age and gender 
and did not differ with respect to ethnicity (for details see Table 1). 
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NPSLE patients showed significantly liigiier disease activity as deter- 
mined by the SLEDAI score (NPSLE: mean = 10.8 ± 5.5, non-NPSLE: 
2.0 ± 2.2, p < 0.0001). There were no significant differences between 
age of disease onset (p > 0.99), disease duration (p = 0.95), and SLICC 
scores (p = 0.18) between the two SLE groups. Mini-mental scores 
were within normal ranges in both SLE groups with no significant differ- 
ences (p = 0.60). For details see Table 1. 

3.2, Conventional MRl analysis 

On the conventional MRl a mild lesion burden was seen in 36.8% 
(7/19), a moderate in 21.1% (4/19), and a severe in 10.5% (2/19) in 
the NPSLE patients. A mild lesion burden was seen in 42.1% (8/19) 
of the non-NPSLE patients, and a moderate lesion burden was seen 
in 15.8% (3/19) of the patients. In the HC group 33.3% of study sub- 
jects (6/18) had a mild lesion burden, 1 subject (5.6%) had a moder- 
ate and 1 subject (5.6%) had a severe lesion burden. There were no 
significant differences in lesion load between the three groups: 
non-NPSLE vs. HC, p = 0.55 (Fisher's exact test); NPSLE vs. HC: 
p = 0.43, non-NPSLE vs. NPSLE: p = 0.61. 

Mild atrophy was present in 8 of the 19 NPSLE patients and in 5 of 
the 19 non-NPSLE patients. A few incidental findings were found 
such as a small pineal cyst present in one NPSLE patient and a small 
meningioma in another NPSLE patient, while a small venous angio- 
ma was present in one non-NPSLE. A few incidental findings were 
also found in the HCs; 4 venous angiomas and one right carotid ar- 
tery aneurysm. 

3.3, Voxel based analysis with threshold free cluster enhancement 

The results of our statistical analyses are summarized in 
Tables 2, 3 and 4, as well as Figs. 1 and 2. Analysis with TFCE re- 
vealed several regions of decreased FA in prefrontal white matter 
in SLE patients, both with and without neuropsychiatric symptoms. 
In total in NPSLE patients 7429 skeleton voxels (at the 1 mm^ res- 
olution used by TBSS) showed significantly decreased FA values 
as compared to HCs. Non-NPSLE patients had decreased FA values 
in 3520 voxels. However, when we directly compared NPSLE and 
non-NPSLE patients there were no significant differences between 
the two groups. 

3.4, Cluster based analysis 

Statistical analyses with cluster forming threshold t > 3.0 revealed 
two regions, with 123 and 110 voxels in the prefrontal white matter 
in each hemisphere, of decreased FA in NPSLE patients as compared to 
HCs (p = 0.01 and p = 0.012). No significant clusters were found 
when comparing non-NPSLE patients to HCs or to NPSLE patients. 
Thus, we obsereed a greatiy reduced sensitivity of cluster-based analysis 



Table 3 

Differences in regional fractional anisotropy (cluster corrected with voxel threshold t > 3.0). 
Comparison of groups Cluster size Coordinates p value 

X y z 

NPSLE < HC 

123 -24 31 5 0.010 

110 24 26 2 0.012 

Non-NPSLE < HC 
NPSLE < non-NPSLE 



HC = healthy control, NPSLE = neuropsychiatric systemic lupus eiythematosus; non- 
NPSLE = systemic lupus erythematosus without neuropsychiatric symptoms. In a first 
step voxels were thresholded at t> 3.0; permutation-based cluster statistics were then 
computed to obtain a multiple comparison corrected p values for each cluster. Coordinates 
refer to voxels with largest values within a cluster. 



Table 4 

Differences in regional fractional anisoti'opy (cluster corrected with voxel threshold t> 2.0). 
Comparison of groups Cluster size Coordinates p value 

X y z 

NPSLE < HC 

1683 24 26 2 0.001 

1382 -24 31 5 0.002 

Non-NPSLE < HC 

536 -6 26 -1 0.027 

NPSLE < non-NPSLE 



HC = healthy control, NPSLE = neuropsychiatric systemic lupus erythematosus; non- 
NPSLE = systemic lupus erythematosus without neuropsychiatric symptoms. In a first 
step voxels were thresholded at t > 2.0; permutation-based cluster statistics were then 
computed to obtain a multiple comparison corrected p values for each cluster Coordinates 
refer to voxels with largest t values within a cluster 



as compared to TFCE. We hypothesized that this might be caused by the 
fact that the large clusters found by TFCE are broken down into smaller 
pieces by the arbitrary cluster-fonning threshold, and that those indi- 
vidual pieces are no longer significant. We then repeated cluster- 
based analysis with a reduced threshold of t > 2.0, which indicated spa- 
tially much more extended regions of significant FA change in NPSLE, as 
coinpared to HCs (two clusters with 1683 voxel, p = 0.001, and 1382 
voxels, p = 0.002), and re-introduced one cluster with significantly 
lower FA values (536 voxels, p = 0.027) in non-NPSLEs as compared 
HCs. For more details see Tables 2 and 3 and Figs. 1 and 2. 

3.5. Multivariate analysis 

Classification analyses revealed a better than chance classification 
when testing all SLE patients vs. HCs (79%, p = 0.001 ), NPSLE patients 
vs. HCs (70%, p = 0.026), and non-NPSLE patients vs. HCs (70%, p = 
0.020). We employed a permutation-based test, which helps to guard 
against accidental double dipping and corrects for imbalanced class 
sizes (e.g. when classifying all SLEs vs. HCs, an accuracy significantly 
higher than 67% was tested). Using the same method, it was not possible 
to classify NPSLE vs. non-NPSLE patients (34% accuracy). 

4. Discussion 

In this study we applied DTI, using both univariate and multivariate 
analysis methods, to determine differences in white matter integrity be- 
tween SLE patients, with and without neuropsychiatric symptoins, and 
a group of HCs. Applying voxel based statistics and using the TFCE 
inethod to correct for multiple comparisons we found decreased FA 
predominantly in prefrontal white inatter in both NPSLE and non- 
NPSLE patients as compared to HCs; we found no significant differences 
between the two SLE groups. However as compared to the HC group 
NPSLE patients displayed about twice as many voxels with significandy 
reduced FA values as the non-NPSLE group. The implications of 
our study are twofold: non-NPSLE patients do show alteration in 
white matter integrity, although to a lesser degree than NPSLE patients. 
Secondly for the categorical decision whether non-NPSLE patients are 
more like NPSLE patients or more like HCs the choice of the statistical 
test, specifically the correction for multiple comparisons plays a critical 
role. 

Our data are in line with other studies describing diminished white 
matter integrity in NPSLE patients (Emmer et al., 2010; Jung et al., 
2010b; Zivadinov et al., 2013), suggesting that directed diffusion is 
reduced in this patient group, which can possibly be attributed to local 
pathologies or a combination of these. Possible mechanisms that have 
been suggested are damage to the myelin, i.e. local demyelination, 
damage to gray matter resulting in secondary axonal loss in the 
white matter, or secondary involvement of the white matter due to re- 
peated small vessel vasculitis and acute vessel wall inflammation 
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(Bosma et al., 2000). Such inflammation of tfie small vessels of the brain 
might then lead to vasculopathy and cause focal hypoperfusion or 
microinfarcts that cause subde changes, not visible on conventional MRl, 
but detectable by DWl or DTI (Hughes et al., 2007; Luyendijk et al., 201 1 ). 

Furthermore our data suggest that non-NPSLE patients, although 
clinically unaffected with respect to neuropsychiatric symptoms, also 
display alterations in white matter structure, in terms of decreased FA 
values in the prefrontal cortex, however to a lesser degree than NPSLE 
patients. Non-NPSLE patients displayed about half as many voxels of 
significantiy reduced FA as NPSLE patients. 

Emmer et al. were the first to apply DTI and TBSS in a cohort of 1 2 
SLE patients and 28 HCs, reporting decreased FA values in various 



regions such as the inferior fronto-occipital fasciculus, the fasciculus 
uncinatus and the posterior limb of the internal capsule (corticospinal 
tract), as well as in the anterior thalamic radiation (Emmer et al., 
2010). Importantly only 7 of these patients fulfilled the ACR criteria 
for NPSLE, so that it is likely that non-NPSLE patients also contributed 
to the group differences. Investigating a study sample comparable to 
ours with respect to sample size and mean age, and applying the same 
methodological approach, with a cluster forming threshold of t > 3.0 
Jung et al. could identify regions of decreased FA only in NPSLE patients, 
but not in non-NPSLE patients (Jung et al., 2010b). Furthermore two re- 
gions showed decreased FA values in NPSLE patients as compared to 
non-NPSLE patients, i.e. the left superior longitudinal fasciculus and 




TFCE 



Cluster- 
based 
t>3.0 



Cluster- 
based 
t>2.0 



Fig. 2. Non-NPSLE patients compared to healthy controls. A significant reduction in regional fractional anisotropy in prefrontal white matter was found in non-NPSLE patients (red voxels). 
No significant differences were found with cluster threshold t > 3.0. 
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the corpus callosum. Importantly, when we applied the same correction 
(for multiple comparisons) method as in (Jung et al., 2010b) our results 
were similar to those reported by Jung et al., i.e. we saw no differences 
between non-NPSLE patients and HCs. For the following reasons, we be- 
lieve that the more profound differences found by TFCE are genuine, 
rather than false positives; first, we could achieve similar results when 
reducing the cluster-forming threshold. We note that there is no sys- 
tematic or "correct" way of setting this threshold, and that it is one of 
the main advantages of TFCE that it does not require such a threshold. 
Its increased sensitivity compared to competing methods is one of the 
benefits claimed by TFCE (Smith and Nichols, 2009); this claim is further 
substantiated by the fact that multivariate analyses, using SVMs and the 
leave-one-out method, were able to classify the two classes (non-NPSLE 
and HC) with better-than-chance accuracy. Given the rather diffuse pat- 
terns of altered FA in SLE, regardless of neuropsychiatric symptoms, we 
did not expect a classifier based on localized features to achieve highly 
accurate classification results. Rather, our motivation to use features 
closely related to the quantities used in statistical testing was to reaffirm 
the significant differences found by TFCE. Indeed, the accuracy to differ- 
entiate between the two SLE groups and HCs was (as expected) modest, 
but still significantiy better than random chance, ( non-NPSLE patients — 
HCs: 70%; and NPSLE patients - HCs: 70%). Differentiation between the 
two SLE groups was not possible (NPSLE patients - non-NPSLE patients: 
34%). Against this background, we tend to think that the results revealed 
by the TFCE method are valid and that reductions in white matter integ- 
rity can be observed in both NPSLE patients and in non-NPSLE patients. 
However, apart from the thresholding procedure other possible differ- 
ences between previous studies and the current might explain, why so 
far no significant differences in the non-NPSLE have been reported; 
these include local patient characteristics and differences in scanning 
parameters. While we were not able to identify features that made 
our sample somewhat special, the current study is the first to apply 
the TBSS algorithm on DTI data that were collected on a 3 T scanner, 
with 15 gradient directions and thus might have caught subtle differ- 
ences that escaped other scanning protocols. Importantly the lesion 
load seen in this study is comparable to that seen in previous studies 
showing a higher percentage of patients with a severe lesion load in 
the NPSLE group, suggesting that both SLE groups are representative. 

Different pathogenic pathways have been proposed for neuropsy- 
chiatric manifestations in SLE and current research tries to identify 
mechanisms that make SLE patients prone to CNS involvement. Often 
neuropsychiatric symptoms occur early in the disease (NPSLE), which 
for some authors speaks for an immune-mediated pathogenesis 
(Steup-Beekman et al., 2013). As such CNS involvement might results 
from a specific antibody profile (Emmer et al., 2010; Matus et al., 
2007; DeGiorgio et al., 2001 ), possibly associated with a leakage of the 
blood-brain-barrier (Huerta et al., 2006) or it might more generally re- 
flect a higher disease activity as compared to non-NPSLE. Given that 
with respect to CNS involvement HCs and NPSLE patients are two 
poles of one spectrum, ranging from healthy to severely ill, our data sug- 
gest that non-NPSLE patients fall somewhere in between. In support of 
this notion white matter lesions, as determined by conventional neuro- 
radiological assessment have been described to occur in up to 60% of 
non-NPSLE patients (Cotton et al., 2004) and cognitive dysfunction 
has been reported to be present in 30% of non-NPSLE patients 
(Monastero et al., 2001 ; Jung et aL, 2012), both strongly suggesting ce- 
rebral involvement despite the absence of neuropsychiatric pathologies 
fulfilling ACR case definitions. One interesting question is how many of 
the non-NPSLE patients in our cohort would have shown cognitive dys- 
function as determined by standardized neuropsychological testing. 
This issue was not fully addressed, and can be seen as a limitation; as 
such the Mini-Mental questionnaire with normal scores in both 
SLE groups is not sensitive to detect subtle, but clinically relevant 
dysfunctions. 

The pattern seen in this study in both NPSLE and non-NPSLE very 
much resembled each other; the other two studies applying DTI and 



TBSS however found partially overlapping, partially different regions 
of diminished white matter integrity. All three studies could identify 
prefrontal areas and parts of the internal capsula, on the other hand 
none of the studies found FA decreases in the occipital lobe. Overall it 
appears that SLE is associated with a rather diffuse pattern of decreased 
FA. Against this background it will be of importance to further DTI stud- 
ies to subclassify the study samples (i.e. NPSLE patients only with one 
symptom, e.g. psychosis, acute confusional state) in order to identify 
more symptom specific FA correlates. 

There are a couple of limitations that need to be pointed out. 
Although great care was taken to come up with study samples, that 
did not differ in group size, age, and ethnicity, and despite SLE 
group performance in the MMSE within normal ranges, and despite 
no significant differences in white matter lesion burden, we cannot 
rule out that some of the differences seen in FA maps are associated 
with differences in behavioral or biological parameters that had not 
been assessed within the study. The characterization and work up 
of patients and HCs in future studies will need to include more pa- 
rameters (e.g. neuropsychological performance, mood) to further as- 
sess the neurobiological correlates of FA alterations. Secondly small 
white matter lesions were found in 8 out of 18 HCs (44%). Small 
white matter lesions have been reported in up to 45% of healthy peo- 
ple (Fazekas et al., 1988; Schmidt et al., 1992), possibly associated 
with a higher risk of experiencing a stroke later in life. As such it 
might be argued that the control group was not ideal. However, 
HCs were healthy at the time point of scanning with no history/evi- 
dence of cardiovascular risk factors or specific white matter disease 
(e.g. leucodystrophies). Including only HCs, that had been deter- 
mined post hoc to have absolutely no lesions, could have induced a 
bias, possibly going along with false positive FA differences between 
groups not related to the SLE status. Importantly white matter le- 
sions and clusters of FA difference were not systematically found in 
the same regions, supporting the notion that small white matter le- 
sions occur independently from SLE associated decreases in FA 
values. 

5. Summary 

Our data demonstrate changes in white matter integrity in terms 
of decreased FA values, predominantly in the prefrontal cortex sug- 
gesting cerebral involvement in NPSLE patients and to a lesser de- 
gree in non-NPSLE patients. These findings go beyond the 
unspecific white matter lesions reported to be present in up to 60% 
of these patients group. Whether this pattern of diminished FA is 
specific to SLE remains to be determined. Clinical and radiological 
work up of SLE patients should be extended to include both neuro- 
psychological testing and DTI in both NPSLE and non-NPSLE patients, 
with the aim to specifically investigate neurocognitive dysfunction 
on the one hand and reductions in white matter integrity on the 
other hand. DTI with a scanning time between 8 and 1 2 min, together 
with advanced analysis methods, is generally suited to become part 
of the routine neuroradiological work up facilitating the investiga- 
tion and assessment of white matter in SLE patients, possibly even 
providing at tool to support or discourage the diagnosis SLE and to 
monitor subtle changes in white matter integrity over time not as- 
sessable through conventional neuroimaging. 
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